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ABSTRACT 
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In the binary radio pulsar system J0737— 3039, the faster pulsar A is eclipsed once 
per orbit. A clear modulation of these eclipses at the 2.77 s period of pulsar B has 
recently been discovered. We construct a simple geometric model which successfully 
reproduces the eclipse light curves, based on the idea that the radio pulses are attenuated 
• by synchrotron absorption on the closed magnetic field lines of pulsar B. The model 

explains most of the properties of the eclipse: its asymmetric form, the nearly frequency- 
^vq ■ independent duration, and the modulation of the brightness of pulsar A at both once 

and twice the rotation frequency of pulsar B in different parts of the eclipse. This 
[ detailed agreement confirms the dipolar structure of the star's poloidal magnetic field. 

"^h \ The inferred parameters are: inclination angle between the line of sight and orbital 

Qlh| plane normal ~ 90.5°; inclination of pulsar B rotation axis to the orbital plane normal 

O ' ~ 60°; and angle between the rotation axis and magnetic moment ~ 75°. The model 

■ makes clear predictions for the degree of linear polarization of the transmitted radiation. 

The weak frequency dependence of the eclipse duration implies that the absorbing 
£> . plasma is relativistic, with a density much larger than the corotation charge density. 

Such hot, dense plasma can be effectively stored in the outer magnetosphere, where 
cyclotron cooling is slow. The gradual loss of particles inward through the cooling radius 
is compensated by an upward flux driven by a fluctuating component of the current, and 
by the pumping of magnetic helicity on the closed field lines. The trapped particles are 
heated to relativistic energies by the damping of magnetospheric turbulence and, at a 
slower rate, by the absorption of the radio emission of the companion pulsar. A heating 
mechanism is outlined which combines electrostatic acceleration along the magnetic 
field with the emission and absorption of wiggler radiation by charged particle bunches. 



1. Introduction 

The double pulsar system PSR J0737— 3039A/B contains a recycled 22.7 ms pulsar (A) in a 2.4 
hr orbit around a 2.77 s pulsar (B) (Burgay et al. 2003; Lyne et al. 2004). Pulsar A is eclipsed once 
per orbit, for a duration of ~ 30 s centered around superior conjunction. The width of the eclipse 
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is a weak function of the observing frequency (Kaspi et al. 2004). Recently, the pulsar A radio 
flux has been found to be modulated by the rotation of pulsar B during the eclipse (McLaughlin 
et al. 2004). The eclipse is longer when the magnetic axis of pulsar B is approximately aligned 
with the line of sight (assuming that radio pulses are generated near the magnetic axis of B). In 
addition, there are narrow, transparent windows in which the flux from pulsar A rises nearly to the 
unabsorbed level. These spikes in the radio flux are tied to the rotation of pulsar B, and provide 
key constraints on the geometry of the absorbing plasma. 

The physical width of the region which causes this periodic modulation is comparable to, 
or smaller, than the estimated diameter of the magnetosphere of pulsar B. Combined with the 
rotational modulation, this provides a strong hint that the absorption is occurring within the 
magnetosphere of pulsar B, which retains enough plasma to effect an eclipse. An alternative 
possibility, which we disfavor, is that the eclipse is caused by absorption in the relativistic wind 
of pulsar A after it is shocked and decelerated by its interaction with the magnetosphere of B 
(Lyutikov 2004; Arons et al. 2004). 

Distinguishing between these two models is greatly aided by a quantitative model of the eclipse 
light curve. In this paper, we construct such a model and show that the light curve is consistent, 
in considerable detail, with synchrotron absorption in the pulsar B magnetosphere. The agreement 
is detailed enough to constitute direct evidence for the presence of a dipolar magnetic field around 
pulsar B. We assume that the absorbing plasma is confined within a set of poloidal field lines 
that is rotationally symmetric about the magnetic dipole axis. Considerations of the sourcing and 
heating of this plasma suggest that its density will exceed the Goldreich-Julian density by a few 
orders of magnitude, and that the absorbing particles will be relativistic. Similar conclusions has 
been reached independently by Rafikov k, Goldreich (2004) under different assumptions about the 
sourcing and heating mechanisms. Note that if the plasma has a finite transverse temperature, 
then particles will be reflected from regions of high magnetic field and trapped on the closed field 
lines. Since the trapped particles cool slowly, a large equilibrium density can be established. 

From this perspective, it is clear that the eclipse model of Lyutikov (2004) and Arons et al. 
(2004) had to make three fine-tuning assumptions. First the wind of pulsar A must be extremely 
dense: about ~ 10 6 times the Goldreich-Julian density at the pulsar speed of light cylinder, in 
contrast with models of pair creation which suggest values closer to ~ 10 3 (e.g. Harding et al. 
2002; Hibschman & Arons 2001). Second, the pulsar A wind must be slow, with typical Lorentz 
factor 7a ~ 10, while at the distance of ~ 1000 light cylinder radii of pulsar A it is expected that 
7A ~ 1000 (Michel 1969). Third, the orbital inclination angle is < 86°; new data indicate that the 
inclination angle is very close to 90° (Cole et al. 2004; Ransom et al. 2004). The model does 
not offer a simple explanation for the observation transparent windows in the eclipse light curves: 
while the position of the magnetosheath is expected to vary by ~ 25% between different phases, 
it is difficult to see how the the sheath can become fully transparent. In addition, the sheath is 
expected to be wider when the magnetic axis is perpendicular to the line of sight while, in contrast 
to this, the eclipse is narrow at such points (McLaughlin et al. 2004). 
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The presence of a such a dense plasma within the magnetosphere of pulsar B is contrary to what 
is usually claimed for isolated radio pulsars. Thus, a source of charged particles — of both signs — 
is required on the most extended closed magnetic field lines. We describe two related mechanisms 
by which the absorbing particles can be supplied by the star itself, or by pair cascades close to its 
surface. 1 First, torsional deformations of the pulsar B magnetosphere, driven by turbulence in the 
surrounding magnetosheath, will generate a fluctuating component of the current, which modulates 
the effective opening angle of the field lines close to the magnetic poles, and the Goldreich-Julian 
current. A second, more subtle, effect involves the pumping of magnetic helicity from the outer 
to the inner magnetosphere. We outline a simple energetic reason why this will occur when the 
outgoing and return currents are not precisely balanced in each magnetic hemisphere. The dipolar 
magnetic field then supports a static twist, and a persistent current flows along magnetic field lines 
which close within the magnetosphere. 

The plasma on closed field lines is quickly heated to relativistic energies. The outer magneto- 
sphere of pulsar B is heated both by the damping of relativistic turbulence, and by the resonant 
absorption of the radio waves from pulsar A. (The resonant absorption of high energy emission 
from pulsar A can occur only very close to the neutron star, where the geometric cross section is 
small.) Since the energy flux in the wind of pulsar A exceeds the energy flux in its radio emission 
by several orders of magnitude, we argue that turbulent heating is likely to dominate - especially 
in regions of high synchrotron optical depth. 

Turbulence in the surrounding magnetosheath generates Alfven waves on closed field lines. 
The non-linear interaction between the waves creates higher-wavenumber modes and, eventually, 
the dissipation of turbulent energy on small scales. Inside the magnetopause, where the pressure 
of the background magnetic field exceeds the particle pressure, the inner scale of the turbulent 
spectrum occurs when the current fluctuations become too large to be supported by the available 
free charges: the plasma becomes charge-starved at small scales (Thompson & Blaes 1998). As a 
result, charges are accelerated electrostatically along the background magnetic field. Bunches of 
such relativistic charges will generate high-brightness radio emission (wiggler radiation) when they 
scatter off the highest frequency Alfven modes. We show that this radiation can self-consistently 
excite the transverse motion of the particles. The conclusion is that self-absorbed radio photons 
that are generated within the optically thick plasma are the dominant heating mechanism. 

The plan of the paper is as follows. In §2 we review the basic parameters of the PSR 
J0737— 3039 system. In the following section, we describe how electric currents may be excited 
in the magnetosphere of a neutron stars, and how the particles which support the current are 
supplied. In §4 we describe particle heating at large synchrotron optical depth, and calculate 
the distribution of synchrotron-absorbing particles that is used to model the eclipse. Section 5 
presents our calculations of the eclipse light curves. In §6 we discuss modification of the model 



x It is unlikely that closed field lines are seeded by particles from pulsar A wind without, at the same time, producing 
unwanted absorption in the much wider region in between the magnetopause and the wind shock (Lyutikov 2004) . 
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due to presence of ions in magnetosphere. The polarization of the radio waves transmitted through 
the absorbing magnetosphere is calculated in §7. In §8 we discuss predictions and compare with 
alternative eclipse models. In the final section 9 we summarize the broader implications of our 
work. 



2. Basic Parameters 

The semi-major axis of the J0737— 3039 system is a ~ 8.8 x 10 10 cm, the orbital period is 2.45 
hr, and the orbital eccentricity is small, e ~ 0.088. At the time of eclipse, the separation of the 
two pulsars is Dab — 8 x 10 10 cm, their relative orbital velocity is V± = 680 km s" 1 transverse to 
the line of sight, and one orbital degree near eclipse corresponds to a time interval of 20.5 seconds. 
The observed eclipses cover A(p ~ 0.035 radians (2 degrees) of orbital phase. The corresponding 
half-width of the eclipsing material, transverse to the line of sight, is (^A(/))Dab = 1-4 x 10 9 cm. 
This is much smaller than the pulsar B light cylinder, which has a radius c/Hg = 1.3 x 10 10 cm. 
The orbital inclination is close to 90° (Cole et al. 2004; Ransom et al. 2004). 

The spin periods and period derivatives, Pa,b and Pa,b, have been measured for the two 
pulsars. Normalizing the moments of inertia Ia,b to 10 45 g cm 2 , the spindown luminosities are 
L A = IaQa&a = 5.8 x 10 33 Ia,45 ergs s _1 and L B = Ib^b&b = 1-6 x 10 30 Ib,45 ergs s -1 (Burgay 
et al. 2003). 

The magnetosphere of pulsar B is truncated, compared with the magnetosphere of an isolated 
pulsar of the same spin, by the relativistic wind flowing outward from pulsar A. The spindown torque 
of pulsar B is therefore modified. The surface magnetic field can be estimated self-consistently only 
by modeling the interaction between wind and magnetosphere. Different effects can contribute 
to the actual torque (Lyutikov 2004; Arons et al. 2004). For example, the open-field current is 
increased compared with that of an isolated dipole, which implies 

L,=AA^M s c(^) 2 (|^) 2 . (1) 

Here I?ns is the average surface dipole field, and Mb is a parameter that rescales the torque acting 
on pulsar B. 

The magnetospheric radius i? mag is determined by the pressure balance between the wind of 
pulsar A and the magnetic pressure of pulsar B evaluated at, 
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Thus, magnetic field at i? mag depends only on parameters of pulsar A wind. Using our fiducial 
numbers we find 
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On the other hand, magnetospheric radius i? mag and surface magnetic field -Bns depend on details 
of the wind-magnetosphere interaction. Using Eqs. (1-2) we find 
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Below we will use these parameters as our fiducial numbers. 
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3. Injection of Plasma on Closed Field Lines 

To calculate the eclipse light curve, one needs the basic plasma properties in the outer magne- 
tosphere of pulsar B - most importantly, the equilibrium density, particle energy distribution, and 
composition. In this and the next section, we motivate our basic assumption that this confined 
plasma has a large synchrotron optical depth and is relativistically hot. The implications of a 
significant ion component of the plasma are addressed in §6. 

One advantage of positioning the eclipsing plasma on the closed magnetic field lines of pulsar 
B is that a large particle density can build up slowly, over many rotation periods. Heating of the 
particles in the outer magnetosphere, where cooling is long, allows them to be trapped through 
the effect of magnetic bottling. Qualitatively, if particles are leaking out of the outer parts of 
magnetosphere on timescale A times the pulsar period, then a source of particles that generates 
the Goldreich-Julian density each period would result in an equilibrium multiplicity A. The typical 
life times of particles in the outer magnetosphere are indeed long, ~ 10 5 spin periods, see eq. (9). 
Since the required density at the edge of the magnetosphere is ~ 10 5 times the Goldreich-Julian 
density (see §5.3), plasma should be sourced at a rate that is comparable with the Goldreich-Julian 
density per period, or higher. This sourcing can be driven either by torsional Alfven waves (AC 
sourcing); or by a pumping of magnetic helicity from the outer to the inner closed field lines (DC 
sourcing) . 



3.1. Electrodynamics of a Perturbed Magnetosphere 

The magnetopause of pulsar B is expected to be strongly turbulent, as particle-in-cell simula- 
tions demonstrate (Arons, Spitkovsky, & Backer 2004, in preparation). Turbulence in the sheath 
will launch (torsional) Alfven waves on the closed magnetic field lines which carry charge and cur- 
rent. The sign of the current alternates on a timescale comparable to the spin period P, which is a 
few times longer than ~ R nmg /c. It will also alternate on a timescale comparable to the flow time 
behind the wind shock (which is comparable to or somewhat larger than P). The depth to which 
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turbulence may be excited in the magnetosphere of pulsar B will depend on more subtle details, 
such as the resonant interaction between compressive and torsional modes. We will not examine 
them in this paper. 

Torsional Alfven waves, launched by turbulence in the magnetosheath, impose a twist on closed 
field lines. If the typical fluctuating component of the magnetic field in the wave at radius -R mag is 
= A(f)Bp, then there is an associated current J ~ cB^ /47ri? mag = A<f>cBp j '47T-R mag . (Here B^ 
is the toroidal magnetic field, and Acj) is the angle through which the poloidal magnetic field Bp 
is twisted.) As the twist in the magnetic field (the Alfven wave) propagates along the closed field 
line toward the star, the associated current density increases. The ratio of the current density to 
the local corotation current density remains approximately constant, 

J ^ A(f> {neR^y 1 _ l ljV i/4 /A^WcosJ^y 1 ^ 



PGJC COS0 B V c / V - 1 / V °- 3 

Here „ 

n-B n B B 

/>GJ = -^ = -cos^— (6) 

is the corotation charge density (Goldreich & Julian 1969), and Op is angle between the local 
magnetic field and fl. The net current is 

r A(f) _ Rla 

I~^ c B m -^. (7) 

z -"-mag 

In order to supply this current, a minimum particle density is required. If at some point the 
actual particle density is below the minimum value |V x B\/4ire, then charges will be supplied 
from the star or created locally through the reaction 2 7 — > e + + e~ . A minimum charge density 
is always present on closed field lines, the corotation charge density (6). When pulsar B is nearly 
an orthogonal rotator (as is implied by our eclipse model; see §5), 9p is close to 90 degrees, and 
the corotation charge density can fail to supply the fluctuating current even for a small twist angle 
A<f> ~ 0.1. In this case, the current can not be supplied by an upward or downward drift of the 
corotation charge; instead charges of both signs will be supplied to the outer magnetosphere. The 
sign of the instantaneous charge flow will fluctuate in tandem with the sign of the current. 

Electrons and positrons cool slowly by cyclo/synchrotron emission in the outer parts of the 
magnetosphere. Beyond the cooling radius, 

1 2/5 



i? coo i = 1.1 x 10 8 



(B NS R^ S ) B 
5 x 10 29 G cm 3 



cm, (8) 



the cyclotron timescale t cyc = Sm^c 5 /2e 4 B 2 becomes shorter than ~ r/c. Only those particles 
reaching the cooling radius can precipitate to the surface. If the particles in the outer magnetosphere 



2 Because pulsar B is older than ~ 10 8 yrs, and its spindown power is small, the density of target X-ray photons 
is probably too small to allow pair creation via the channel 7 + 7 — > e + + e~ 
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are isotropized, then the cooling fraction is small (§4.3), 
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The upward flux of particles is shut off when the local density at -R coo i exceeds the density \J\/ec = 
AjraGj(-Rcooi) that is needed to supply the entire current close to the neutron star. On this basis, 
we obtain an estimate of the equilibrium multiplicity of electrons (or positrons) 
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(10) 



in the slow-cooling region. Here r is the distance from pulsar B at the magnetic equator. Notice 
that this expression is independent of the electron temperature in the outer magnetosphere. It also 
implies a characteristic number density, 
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that is approximately independent of r > -R coo i- 



An independent upper bound on the multiplicity comes from the requirement that the particle 
pressure not exceed the magnetic dipole pressure: 
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Electrons/positrons with a density A cq have, equivalently, a maximum temperature 
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The supply of charges to the magnetosphere can be greatly augmented by multiple stages of 
pair creation (e.g. Hibschman & Arons 2001). This requires the formation of large gaps (regions 
of uncompensated parallel electric field) in the inner magnetosphere. From the perspective of the 
near-field electrodynamics of the neutron star, the twist in the magnetic field is modulated only 
very slowly (over ~ 10 4 times the light-crossing time of the star). This part of the circuit will 
therefore bear some resemblance to the polar cap of an ordinary, isolated radio pulsar. 

Electric gaps can appear higher in the magnetosphere, if the particles at a given point have 
a fixed sign and a charge density equal to />gj- These gaps will form along the surfaces where 
ft ■ B = (e.g. Holloway 1977). The magnitude of the electrostatic potential drop that builds 
up across this gap during one torsional oscillation would approach A<I> ~ eBRAcj). However, the 
presence of a dense, polarizable plasma would largely prevent these gaps from forming on the closed 
magnetic field lines, outside the cyclotron cooling radius (8). The gaps would also be suppressed 
if the magnetosphere supported a persistent twist with a magnitude greater than ~ Jlr/c, a 
possibility that we entertain in Appendix A. 
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3.2. Comparison: Absorption of Particles from the Wind of pulsar A 

It is also worth considering the wind of pulsar A as a competing source of absorbing particles 
in the magnetosphere of pulsar B. Transfer of particles from wind to magnetosphere will occur 
rapidly through reconnection between the oscillating magnetic field that is advected by the wind, 
and the forward surface of the magnetosphere (e.g. Thompson et al. 1994). The particle density 
in the pulsar A wind at the position of pulsar B is 



Xa Ia^a^a . in _2, -3 ( -* a\ 

^ = i^Vi^ = 4x10 A ^ cm ■ (14) 



Here Xa is the multiplicity of pairs created by a cascade on the open field lines of pulsar A. 
Comparing with the characteristic charge density ng = A eq nGj,_B = A eq Jl.B-B(-R ma g)/27rec in the 
magnetosphere of pulsar B, at a distance -R mag from the star, one has 



^ ^(-Rmag) £ A C q / Rmag 



5 



ua Xa \4 x 10 9 cm 



5 x 10 29 G cm 3 



(15) 



The density of particles trapped from the wind of pulsar A will not, generally, exceed the density 
in the wind itself: otherwise the return of particles from the magnetosphere back to the wind will 
balance the gain. One expects that A cq > 10 5 > A^, and so we conclude that the dominant source 
of absorbing particles in the magnetosphere of pulsar B is likely to be the particles pulled outward 
from its inner magnetosphere. 



4. Particle Heating in the Outer Magnetosphere 

There are two principal sources of free energy in the outer magnetosphere of pulsar B: the 
relativistic wind emitted by pulsar A, and the radio emission of pulsar A. Both can be effective 
at heating trapped particles to relativistic energies. The wind energy flux is, nonetheless several 
orders of magnitude larger (in spite of the uncertain effects of wind collimation and beaming). 
In addition, the equilibrium plasma density estimated in §3 is too large to allow effective heating 
through radio absorption throughout the bulk of the pulsar B magnetosphere. 

The kinetic power incident on the magnetosphere (radius Rmug] eq. 4) at a separation Dab = 
8 x 10 10 cm is 

Pwi„d * (^) 2 La = 4 x 10 30 A- B 1/2 (/a,45 W /2 ergs s" 1 . (16) 

The radio energy flux incident on pulsar B is vF v = (D J0737 / Dab) 2 ( I/ -^l/)obs — 1-2x10 (.DJ0737/O.6 kpc) 2 
ergs cm -2 s" 1 at 1.4 GHz, where (i ? 1 y) bs ~ 1-6 mJy is the observed flux at this frequency, and 
£>J0737 ~ 0.6 kpc is the estimated distance of the PSR J0737-3039 system Lyne et al. (2004). 
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The incident radio power is smaller by 
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4.1. Heating by External Radio Photons 

In spite of the small total radio power output, the absorption of the radio waves can have an 
important influence on the kinematics of trapped electrons. [See Rafikov & Goldreich (2004) for 
an independent and more detailed analysis of this heating mechanism.] For present purposes, we 
assume the existence of a flux of electrons (and possibly positrons) moving trans-relativistically 
along the closed poloidal magnetic field line, away from the star. We now show that these particles 
will be heated sufficiently to mirror and become trapped in the outer magnetosphere. The trapped 
particles are then further heated by the radio beam. 

A non-relativistic electron absorbs (unpolarized) radio photons with a cross section <7 cyc (^) = 
(ire 2 /2m e c) [1 + cos 2 k\5(v — fB,e)- Here VB,e = eB /2irm e c and k is the angle between the direction 
of magnetic field and the line of sight. The energy absorbed by one particle from the radio beam 
at a radius r (outside the cooling radius (8)) is 



rdE 
v dt 



a cyc (y)F v dv, 



(18) 



when the particle motion is trans-relativistic. This works out to 
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10 27 ergs s" 



keV, 



(19) 



when 1 - Vl-^ 2 /c 2 = 0(1). This energy is absorbed far outside the cooling radius, and so the 
particle mirrors soon after it begins to return to the star. 

Particles trapped in the magnetosphere will continue to be heated. The rate of heating at any 
one position in the magnetosphere depends on the column of intervening particles, and the pump 
spectrum. Using synchrotron absorption coefficient for a thermal distribution of particles at low 
frequencies we find that the time to heat the particles up to a temperature T 3> m e c 2 /&b is short: 
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(20) 



where f pea k,9 = ^peak/10 Hz. Balancing optically thin synchrotron cooling at a rate lQ(k-e,T / m e c ) 
(n e aTc)(B 2 /8ir) per unit volume, with synchrotron heating at a rate a u {vF v ) Upc&]!i per unit volume, 
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one obtains an equilibrium Lorentz factor 
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This is fairly close to the limiting temperature (13) obtained by balancing the particle pressure in 
the outer magnetosphere with the magnetic pressure. 

4.2. Electrostatic Heating and Thermalization of an Optically Thick Plasma 

The wind energy of pulsar A that is incident on the magnetosphere of pulsar B will be con- 
verted, with some efficiency, to internally generated radio waves. Even if the radio output of the 
magnetosphere L mag is much weaker than that of pulsar A, its heating effect can dominate by the 
factor 4(Dab / Rma,g) 2 L mag / La ~ 2000(i? mag /4 x 10 9 cm)~ 2 L mag /LA- So we consider the possibility 
that the magnetospheric plasma is, itself, a source of low frequency photons {y < 100 MHz) which 
are created and absorbed in situ. The unpulsed emission from the PSR J0737— 3039A/B system 
is several times brighter at 1400 MHz than is the combined pulsed emission of the two neutron 
stars (Burgay et al. 2003). It is possible that some of this unpulsed emission is generated in the 
magnetosphere of pulsar B. 

Internal heating has another advantage over the absorption of radio photons from an external 
pulsar: it is more effective when the seed particles move relativistically along the magnetic field. The 
freshly injected particles will only resonantly absorb photons of a very low frequency, u ~ 7|j~ 1 ^B, e - 
At this frequency, the external radiation sees a very large optical depth (from the previously injected 
particles which have already acquired large perpendicular energies) . 

When the density n e of electrons (and positrons) greatly exceeds the corotation charge density, 
torsional waves generated with a wavelength fciir ~ 1 will couple non- linearly to higher-frequency 
waves. If the energy density in the light charges is still small compared with that of the background 
magnetic field, the inner scale of the resulting turbulent spectrum is determined by balancing the 
fluctuating current density with the maximum conduction current that can be supported by the 
plasma (Thompson & Blaes 1998). At high frequencies, the cascade is expected on general grounds 
to be strongly anisotropic. The wavenumber kj_ of a wavepacket perpendicular to the background 
magnetic field will be much greater than the parallel component k\\ , and related to it by k±5B ~ k\\B 
(Goldreich & Sridhar 1997). Balancing k±SB ~ 47rn e e, and relating n e to the corotation charge 
density through n e = X cq (QB /2nec), gives 



Relativistic electrons moving anti-parallel to such an Alfven wavepacket with Lorentz factor 




(22) 



7|| will emit photons of a low frequency (wiggler radiation), 

lo a ~ l\k\\c 



(23) 
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(Fung & Kuijpers 2004). Large parallel Lorentz factors are most easily achieved by particles which 
are freshly injected with small perpendicular energies into the outer magnetosphere, and then 
electrostatically accelerated in regions of high turbulent intensity. These same charges will absorb 
the wiggler photons at their rest-frame cyclotron resonance (lab frame frequency eB /2Trm e cy^) if 



— i— )=—(—). (24) 
2A cq \m e cn B ) A cq V^mag' 



This condition is easily satisfied for the multiplicities (A cq ~ 10 5 ) that we encountered in §3. It is 
straightforward to check that the implied wiggler frequency (23) lies comfortably above the plasma 
frequency at density A cq . It works out to 

^ = 3.»f^V /3 (l^T MH- (25) 
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The radio power needed to heat the freshly injected electrons is a small fraction of the net power 
dissipated in the magnetosphere. For example, if the particles carry a current J that supports a twist 
A(f> in the background magnetic field, then their kinetic power (before electrostatic acceleration) is 

4x 1023 fe) (£;)"' < 26 > 

The minimum radio power that is needed to heat these particles is no larger than L\^ n . (It should, 
nonetheless be emphasized that some bunching of the radiating particles within the turbulent 
plasma is required to exceed this requirement.) 

Relativistic motion of the injected particles reduces the radio power even further. The particles 
start in their lowest Landau states close to the star. As they move out beyond the cooling radius (8), 
they will begin to decelerate as they absorb photons. The Lorentz factor parallel to the magnetic 
field is halved when the energy of the absorbed photons in the particle rest frame is AE ~ 7_|_m e c 2 . 
Here 7j_ = (p±/m e c) 2 + 1. The energy of the absorbed photons in the star's frame is smaller by 
a factor ~ 7jj" 1 . The net density of the absorbed photons is therefore smaller than the beam energy 
density, 

XJ 1 ~ — n beam m e c 2 = -o^bcam oc 7I? 1 - (27) 
7|| 7y 11 

The equilibrium electron temperature can be much higher in this situation, than we found by 
assuming the external radio photons from pulsar A to be the sole radiative pump. Balancing the 
rate of turbulent heating with the incoherent synchrotron output of the thermalized particles, and 
assuming that a fraction eturb of the wind energy density Bmag/^ 71 " ^ s damped on a timescale r/c, 
one finds 

^S) 2 = 2 * iol M^) 8 < 4 < 28 > 
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So, for example, if the energy density in torsional Alfven waves is a fraction ea of the wind energy 
density at a distance r from pulsar B, then the three- wave damping timescale is ~ (6 B / B)~ 2 (r / c) ~ 
s^\B/B mag ) 2 (r/c). One finds 

eturb = A (ir-) = £ A (-^—) ■ (29) 

\ "mag / \ -rtmag / 

The implied equilibrium temperature is 

k * T * - 2S0, A ( ^V"Vr ( TTT ^ = Y (30) 



m e c 2 \W 5 J B V!-5xl0 9 cm 

For ea ~ 0.1 this estimate is consistent with Eq. (13). 



4.3. Density Distribution of the Heated Particles 



Suppose that at the edge of the magnetosphere, located at radius -R mag , relativistic particles 
with a typical Lorentz factor 70 are isotropised. Particles captured from the wind of pulsar A, and 
particles injected into an optically thick plasma from below (§4.2) are both expected to satisfy this 
criterion approximately. Neglecting for the moment the effects of cyclo/synchrotron emission, the 
pitch angle ip can be related to the value VWg at a radius r 



Rmag from the conservation of the first 



adiabatic invariant p\/B and the constancy of the Lorentz factor 70 . [Here pj_ 
is the is particle momentum perpendicular to the magnetic field.] One finds 



sin?/' 



B 



B 



mag 



1/2 



sin tp- 



mag • 



sin^ (7q — 1) l l 2 m e c 



(31) 



Reflection occurs when tp = it/ 2, where the field strength is 

Bred 1 



B 



nia« Sin 2 V'mag 



(32) 



Near the axis of the dipolar field, the reflection radius is i? r efl/-Rmag — sin 2//3 VWg • The particle 
spends only a short time ~ R rc n/c near the reflection point, as may be seen from the solution for 
the speed of the particle parallel to the magnetic field, 



;|| = vq cos ip = vq 



B 



B 



mag 



I 1/2 



sin tp 



mag 



(33) 



Here vq/c 



1 



7o 



The density of particles with initial pitch angle tp ma , g 
particle flux, 



is obtained from conservation of the 



dn v\\ 



dn 



Vq cos ip- 



mag 



dlP 



mag 



B 



dlp : 



mag 



mag 



B 



(34) 
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Integrating over initial pitch angles up to a maximum value ^ max = arcsin(i? ma g/i?) 1//2 gives the 
local density of particles, 



dn 

n 



mag 



mag 



f*™** (g/gmag) COS ^ mag Sin^magd^mag _ 

[l-(S/5 mag )sin 2 ^ ma g] 1 



(The factor of 2 in front of this expression accounts for the reflected flux of particles.) Thus, the 
density of particles is approximately constant, in spite of the strong convergence of the magnetic 
field lines toward smaller radius. 

The adiabatic invariant is no longer conserved if a particle reaches deep enough into magne- 
tosphere that it cools significantly. This occurs for electrons (and positrons) inside the radius (8). 
The consequences for the equilibrium density of particles are discussed in §3.1. 



5. Geometrical Model of Eclipses 



In this section we assume that closed field lines of within the pulsar B magnetosphere are 
populated by relativistically hot plasma. We calculate the synchrotron optical depth over a large 
number of lines of sight, taking into account the three-dimensional structure of the magnetosphere. 
We also assume, for simplicity, that the magnetic field is described by a vacuum dipole, and that 
the absorbing plasma is truncating outside some radius. This simple model can reproduce all the 
salient features of the eclipse light curve, including an essential part of the eclipse phenomenology, 
the weak frequency-dependence of the eclipse duration and dependence of phases of B. The optical 
depth to synchrotron absorption at all frequencies (especially at the highest) has a strong gradient 
near the eclipse boundary, and quickly becomes large (t„ > 1) over a small range of orbital phase. 
This can be achieved if the absorbing particles are relativistic and can absorb in a wide frequency 
range. 



5.1. Synchrotron Absorption 

We assume that the closed field lines are populated by relativistic electrons with a thermal 
distribution at a temperature k-BT e /m e c 2 ~ 10 (see eq. [21]). In this case, the peak of the 
synchrotron emission is at a frequency v pea k ~ 4(fcBT e /m e c 2 ) 2 VB,e- The cyclotron frequency is 
VB,e{Bms.g) ~ 2 x 10 7 Hz at the edge of magnetosphere, and increases inward. Therefore radio 
waves propagate in the low frequency regime, v <C f pC ak) when the observing frequency is in the 
gigahertz range. 

The eigenmodes of the electromagnetic wave are linearly polarized when ions are absent from 
the plasma. The synchrotron absorption coefficients of the two eigenmodes are then (Rybicki &; 
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Lightman 1979) 



a 



(1,2) 



2 d 
de 



n(e) 



The polarization-averaged absorption coefficient is 



1 

a v = - 



a 



(i) 



(2) 



(36) 



(37) 



In these expressions, the spectral power density emitted by a single particle is denoted by 

\^3e 3 B sin k 



p(1.2) 



47rm e c 2 



" / K b/Z {rj)dri ± vK 2/z (v) 



v = 



2v 



m e c 2 J 3(sm k)i/b,< 



VB,e 



2-Kvn e c 



(38) 



and we assume an isotropic distribution of pitch angles. The electromagnetic wave is absorbed on 
particles with pitch angle equal to the angle k between B and the direction k of wave propagation. 
In polarization state (1), the electric vector is orthogonal to the /c-B plane (the + sign in eq. 
[36]); whereas in polarization state (2), it sits in the /c-B plane (the — sign in eq. [36]). Also e is 
the energy of the emitting particle, ra(e) is particle distribution function, and K 5 / 3 and if 2/3 are 
modified Bessel functions. 



For a thermal distribution n(e) oc e e / kBT characterized by a temperature T and total density 
n e , the absorption coefficient below the peak frequency is 



Air 2 en e 



a 



(i) 



3 7 / 3 i? sin k 
4 



a 



(2) 



2 



8111 K 



(39) 



5.2. Eclipse Light Curve 

We introduce a Cartesian system of coordinates x, y, z centered on pulsar B (see Fig. 1). The 
plane of the orbit is taken to coincide with the x — y plane, from which the line of sight is offset 
by a vertical distance zq (eq. (67)). The observer is located at x — ► 00. The spin axis of pulsar B 
is inclined at an angle 6q to the orbital normal, and at angle with respect to the x — z plane. 
The magnetic moment of pulsar B has a magnitude [Ib = (BnsRns)b> an( ^ ^ s inclined at an angle 
X with respect to f2#, To calculate the intensity of the transmitted radiation, we need to work out 
the magnetic polar angle # M at each position x = [x, y(t), zo] along the line of sight. 3 Here 

(40) 




3 For the purposes of an initial calculation, we neglect the possible presence of a toroidal magnetic field (Appendix 

A). 
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is expressed in terms of the semi-major axis a cm and the orbital phase <f>. The distance from a 
given point to pulsar B is 

r(x, t) = ^x 2 + y 2 (t) + z 2 . (41) 
The magnetic polar angle of coordinate r is determined from the unit vector fr(t) parallel to 

cos0 M = ^. (42) 
The components of fr are easy to write down in a coordinate system aligned with H^, 

fr% = smxcos(tt B t); fry = smxs'm(n B t); fr% = cosx- (43) 
Transforming to the observer's coordinate system gives 

frx = {frx cos On + sin On) cos fa - fry sin fa; (44) 
fry = fry cos fa + (fr^cosOn + fr%sm9n) sin fa; fr z = fr% cos 6>q - fr% sin 6q . (45) 
The strength of the magnetic field at position x is then given by 



Vl + 3cos^ 

B = ~ ^3 »B- ( 46 ) 

The condition that a given point is located on closed field lines is that the maximum extension of 
a field line is less than i? mag : 

r 

sin 2 0„ 

We also need to know the angle between the line of sight and the local direction of the magnetic 
field: 

rns K _ B - _ 3cos^(x/r)-^ 

COSK " B ~ (l + 3cos2^)V2 • ^ 

For calculation of polarization properties, the electric vectors of the two eigenmodes are parallel to 
the unit normals 

B y z-B z y B y y + B z z 

1 (B-j + £f)l/2' 2 (£2 + £2)1/2 ^ ^ 

which lie in the plane of the sky. 

We first calculate the total intensity of the transmitted radiation. Our examination of polar- 
ization effects is deferred to §7. The optical depth at time t, frequency v, and impact parameter 
zq is obtained by integrating the polarization-averaged absorption coefficient (39) along the line of 
sight, 

/oo 
a u [x,y(t),z ]dx. (50) 
-oo 

We normalize the density n e of absorbing particles to a characteristic Goldreich-Julian density at 
the magnetospheric boundary (B = £ mag ): 

n e = A mag nGj( J Rmag) = A mag ~T^T = ' 2 Ama § cm ~ 3 - ( 51 ) 
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(This is only a convenient normalization: in fact, pis and Qb are nearly orthogonal in our best 
eclipse model.) We also normalize the line-of-sight coordinate and field strength to R mag and -B mag 
(eq. [4]). For our fiducial parameters, the optical depth is 

, 4.5 x 1(T 6 t( x \ , / k B T e \~ 5/3 /Bsin/e\ 2/3 . , 

Here z^ghz is the wave frequency normalized to 1 GHz and we have set I a = Ib = 10 45 g cm 2 . In this 
expression, the parameter Mb encapsulates the uncertainty in the magnetospheric radius through 
the normalization of the spindown torque acting on pulsar B (eq. [4] ) . The code also tests whether 
the local cyclotron frequency satisfies condition fB,e/v > (7e) = (3A;BT e /m e c 2 ), and whether a 
given point is located outside of the cooling radius (8). Both of these effects are important only for 
extremely small impact parameters zq. 

As a simple prescription for density distribution we assume that only those field lines are 
loaded with absorbing particles for which maximum extension (calculated at the magnetic equator) 
lies within a prescribed range, 

-^abs- < -Rmax < Rabs+- (53) 

Motivated by §4.3, the plasma density n e and temperature T e are assumed to be constant along 
each field line which lies in this range. We further assume that n e does not vary between field lines. 

The onset and termination of the eclipse are then determined by the physical boundary of the 
absorbing plasma. We find that i? a bs+ must be fixed at a value ~ 1.5 x 10 9 cm in order to reproduce 
the width of the observed eclipses. The optical depth must quickly become large on lines of sight 
passing just inside the plasma boundary (at a radius ~ i? a b s +)- The implied multiplicity is large, 
Amag ~ 10 5 , but still consistent with our estimates in eqs. (10) and (12). 

One expects a dipole to be only a rough approximation near the plasma boundary when 
-Rabs+ ~ -Rmag- However, the dipole pressure rises rapidly inward, and so the dipole approximation 
becomes increasingly accurate as -R a bs+/-^mag becomes smaller. Modest deviations between the 
model and the data near the edges of the eclipse could be used to probe the distortion of magnetic 
field lines from a true dipole. 



5.3. Successful Eclipse Parameters 

The calculation of the eclipse light curve requires a choice of several parameters: the angles 
<f>n, X-, electron temperature T e and density multiplicity A mag , inner and outer radii i? a b s ± and 
impact parameter zq . Our choice of the parameters was guided by a qualitative sense of how they 
would influence the shape of the light curve, in the following key respects: 

• the eclipse is asymmetric, with the ingress shallower and longer than the egress; 
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• near ingress, transparent windows appear twice per rotation of pulsar B; 

• near the center of the eclipse, the transparent windows appear once per rotation; 

• near egress, the transparent windows are not well defined; 

• the duration of eclipse at a given rotational phase of pulsar B depends on the phase of B, 
being larger when the magnetic moment is along the line of sight; 

• the duration of eclipse at half maximum flux is ~ 30 seconds, while the full duration is ~ 40 
seconds. 

The modulation of the radio flux during the eclipse is due to the fact that - at some rotational 
phases of pulsar B - the line of sight will only pass through open magnetic field lines where is 
assumed to be negligible. One of the main successes of the model is its ability to reproduce both 
the single and double periodicities of these transparent windows, at appropriate places in the eclipse. 
This requires that /7b be approximately - but not quite - orthogonal to fi#. To understand this 
result geometrically, consider Figs. 2-3. The rotation of /7b out of the plane of the sky brings 
the observer's line of sight to within a small angle of the magnetic pole (in the favored geometry). 
The cross-sectional area of the absorbing plasma is maximized at this moment, and it projects an 
ellipsoidal shape on the sky. On the other hand, when /7b sits in the plane of the sky, the absorbing 
plasma fills a region bounded by a set of closed dipole field lines (Fig. 3). The shifting inclination 
of /7b through one half a rotation of pulsar B allows the line of sight to intersect this dipolar region 
once or twice. The proportions of the eclipse in which either periodicity is observed depend on the 
the impact parameter zq and the extent to which \xb is nearly - but not quite - orthogonal to Qb- 
A major success of the model is that it can reproduce the relative sequence and duration of these 
distinct periodicities. 

Our simulated light curve is displayed in Figs. 4 and 5. Variations in the parameters of the 
model modify the light curve in the following ways. A large asymmetry between ingress and egress 
results from a combination of finite 9q (inclination between the spin axis of pulsar B and the orbital 
plane normal) and finite zq (impact parameter). The asymmetry is largely erased if (pa is in the 
range ±60° (the rotation axis is out of the plane of the sky). There can be one or two transparent 
windows per rotation of pulsar B, depending on zq and x (the angle between rotational angle and 
magnetic moment). If x differs considerably from tt/2, then both eclipse center and egress show 
strong modulation on a single rotational period. A very small value of -R a bs- results in a complete 
eclipse in the center. Since v > vb^ in most or all of the eclipse, the inner boundary corresponding 
to VB,e/lcV > 1 is never reached. 

Based on a sample of many eclipse light curves our best fit parameters are: zq — —7.5 x 10 8 cm, 
9q ~ 60°, 4>n — —90°, x — 75°. The orbital inclination is therefore predicted to be close to 90.55°. 
This is somewhat larger than that inferred by Cole et al. (2004) (90.26° ± .13°), but consistent 
with the estimate of Ransom et al. (2004) (88.7° ± .9°). Note that there is a symmetry in these 
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parameters, zq — > — zo and 0q — > —(frn, which preserves the shape of the light curve. In view of 
results of Cole et al. (2004), we have chosen zq < 0. 

The dependence of the eclipse profile on the parameters z , 9^ and % is illustrated in Fig. 8. 
Overall the best constrained parameter is zq ~ —7.5 x 10 8 cm, which must lie within a range a 
range ~ ±10% of this value. There is however, a degeneracy between 9a and x which produces 
nearly identical light curves if both angles are varied up or down by the same factor. The results 
are not sensitive to \(f> n + 90° | < 30°. 

The required minimum multiplicity is A mag ~ 3 x 10 5 Vq^ z (fcBT e /10 m e c 2 ) 5//3 to reduce the 
transmitted radio flux by > 90% at an orbital phase of ~ 90.5°. This electron density is, in fact, 
close to the estimates (10) and (12). The full duration of eclipse is about 40 sec. The size of the 
eclipsing region is i? a bs+ = 1-5 x 10 9 cm, so that the absorbing plasma must be truncated well 
inside the expected radius of the magnetopause. We comment on the significance of this result in 
the next section. 

The model readily reproduces many fine details of the eclipses. It explains the modulation of 
that is observed at the first and second harmonics of the spin frequency of pulsar B, and the deepen- 
ing of the eclipse after superior conjunction. The average eclipse duration is almost independent of 
frequency when the multiplicity is larger than A mag ~ 3 x 10 5 (Fig. 6). Eclipse is broader when the 
magnetic moment of pulsar B is pointing closest to the observer, just as is observed (McLaughlin 
et al. 2004). Figure 7 shows the eclipse profile averaged over different orientations of magnetic 
moment. When Pb is pointing towards the observer the "doughnut" of the closed field lines is seen 
nearly face-on, producing broad eclipses; when fig is in the plane of the sky the "doughnut" is seen 
edge-on, producing narrower eclipses. 

It should also be noted that the details of eclipse ingress and egress are not fully reproduced 
by the model: in particular the rise in the radio flux at egress is sharper than observed (especially 
if A mag is high enough to give the eclipses a weak frequency dependence). This could be due to 
deviations of the actual magnetic field from our assumption of a pure dipole; or be an artifact of 
our assumption of a sharp plasma boundary. We examine the effects of a smoother plasma profile 
and a variable electron temperature profile in §6. 

5.4. Eclipse Duration 

Our eclipse calculations show that the optical depth of the absorbing plasma undergoes a sharp 
drop at a distance i? a bs+ — 1-5 x 10 9 cm from pulsar B. This is about 2.5 times smaller than the 
expected radius of the magnetopause, -R mag — 4 x 10 10 cm (assuming that the torque parameter 
Mb is unity in eq. (4)). The disagreement with the maximum lateral extension of the closed field 



lines 4 is closer to a factor of 4 (e.g. Lyutikov 2004). 

The synchrotron optical depth (52) reflects both the geometrical distribution of absorbing 
particles, and their temperature profile. One possibility is that i? mag has been overestimated, 
and that -R a bs+ reflects the actual boundary of the pulsar B magnetosphere. Alternatively, the 
suppression in r v could be due to particle loss and overheating in the outer magnetosphere. Let us 
consider these possibilities in turn (we clearly favor the latter). 

Recall that the torque parameter Mb is presently unknown. The surface magnetic field of 
pulsar B could be smaller than is implied by the the simplest estimate of the spindown torque 
{Mb ~ 1; eq. (1)). The torque formula (1) assumes an increase in the opening of the dipolc 
field lines that is expected from simple geometry. But other effects are at work. The pressure is 
distributed asymmetrically about pulsar B, so that even if its spin were in corotation with the orbit, 
it would be still be torqued by the pulsar A wind (the Magnus torque, e.g. Thompson et al. 1994; 
Arons et al. 2004). In addition, a strong surface resistivity at the interface between the pulsar B 
magnetosphere and the shocked wind of pulsar A would result in an enhanced spindown torque, 
due to the resistive dragging of the poloidal field lines of pulsar B. 

There is a clear upper bound on the torque and spindown luminosity of pulsar B: 



which comes from the limit of comparable toroidal and poloidal fields at the magnetospheric bound- 
ary. Eq. (54) implies that Mb < c/(f2#i? ma g) ~ 10. If Mb saturated this bound, the magnetospheric 
radius i? mae would be 



The lateral radius of the closed field lines is 50 percent larger, or 3.6 x 10 9 cm, which is still a 
factor of 2 too large. Thus, one cannot explain small eclipse duration as being entirely due to large 
torque/small magnetic field of pulsar B. 

The short eclipse duration may also reflect the loss of plasma from the outermost closed 
field lines, resulting from a change in the topology of the field lines. Experience with planetary 
magnetospheres suggests that the boundary between open and closed field lines can vary by a 
factor ~ 1.5 in polar angle (a factor ~ 2 in maximum dipole radius) from the windward side to the 
leeward side (e.g. Kabin et al. 2004). Indeed, the relativistic PIC simulations of Arons et al. (in 
preparation) show that the plasma is very dynamic in the vicinity of the separatrix between open 
and closed field lines. 



4 Note that the shock which decelerates the wind of pulsar A will have a yet larger transverse dimension (by up a 
factor of 2 or so (Lyutikov 2004; Arons et al. 2004). This disfavors synchrotron absorption in the shocked wind of 
pulsar A as the explanation for the eclipse. 




3 



(54) 




1/3 



(55) 
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If the magnetic field in the pulsar A wind alternates in sign over the pulsar A period (e.g. the 
wind is 'striped') then reconnection can occur over most of the magnetopause surface. On the other 
hand, if the direction of the wind magnetic field is constant in each rotational hemisphere of pulsar 
A, then reconnection will occur only on particular magnetospheric field lines that happen to run 
counter to the wind field. Trapped particles will, nonetheless, be effectively redistributed over the 
bulk of the magnetosphere by azimuthal drift. The drift timescale for particles of Lorentz factor j e 
is ~ i?^ iag i/B,e(-Bmag)/7eC 2 ~ 2 x 10 5 s, which is shorter than the residency time ~ (^_b/ coo i) 1 ~ 
5 x 10 5 s (eq. [9]). As a result, the loss of particles by reconnection in the magnetotail would reduce 
the particle density on the windward side of the magnetosphere as well. 

The loss of particles from the outer magnetosphere is likely to cause overheating of the re- 
maining particles, which reduces the synchrotron absorption coefficient (39). A similar effect will 
occur if the trapped plasma has a large ion component, since the density (56) of the electron-ion 
plasma is suppressed in the outer magnetosphere. This variant of the eclipse model is examined in 
the following section. 

Two other explanations for a short eclipse duration can be discounted. The angle between 
the line of sight and the orbital plane could be larger than is claimed by Cole et al. (2004) and 
Ransom et al. (2004). Increasing the inclination to ~ 2° would reduce the eclipse duration by a 
factor of two. However, our modeling of the eclipse profile places tight constraints on the orbital 
inclination which disfavor this possibility. 

Another unlikely explanation for the narrow eclipse width is that the moment of inertia of 
pulsar B - but not of pulsar A (eq. [4]) - is much smaller than neutron star models suggest. More 
than an order of magnitude decrease in Ib is required, which is not possible for (the lower mass) 
pulsar B even if it were composed entirely of u-d-s symmetric quark matter. 

6. Eclipses by an Ion-Supported Cloud 

We now examine the case where the absorbing plasma is composed largely of electrons and 
ions. In this variant of the eclipse model, the plasma is more centrally concentrated around pulsar 
B than was assumed in §5; and allowance is made for a higher electron temperature in the outer 
magnetosphere. 

The ions have two properties which influence the plasma distribution. First, they have a much 
stronger gravitational binding to the neutron star surface than do electrons and positrons; and, 
second, they cool much more slowly in the outer magnetosphere. These properties of the ions 
have opposing effects on the equilibrium electron density in the outer magnetosphere. On the one 
hand, the large gravitational binding energy makes it energetically favorable to force some of the 
slow-cooling electrons downward through the cooling radius (with a parallel electric field) if their 
temperature is k^T e < (A/Z)GM]sss,Bfnp/ Rns,b ~ 200 (A/Z) MeV. Here Am p is the ion mass, and 
Ze its charge. 
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On the other hand, if the ions are themselves heated at the same time as the electrons, they 
will prevent the lighter charges cloud from collapsing through the cooling radius. The principal 
heating mechanisms which we examine in §4 involve the absorption of low-frequency transverse 
photons (radio waves). For example, heating of protons by such a mechanism will be suppressed if 
the proton cyclotron frequency ujb, p = eB/m p c sits below the plasma frequency of the (relativistic) 
electrons, = (2-7rn e e 2 c 2 jk^T^I 2 . This gives a characteristic electron density 

_ (kBTA B 2 



n e = j - J- ( 56 ) 

\m p c 2 J 2tt m p c 2 



The corresponding lower bound on the multiplicity of the neutralizing electrons (mass m e ) is 

(Sns^ns)^ 



e .i 180 f k B T e /m e & 
eq cosOb 1 10 



5 x 10 29 G cm 3 



4 x 10 9 cm 



(57) 



This can be competitive with expression (10) somewhat inside the magnetopause, e.g. at r < 1 x 10 9 
cm. 

From the requirement that the electron pressure not exceed the magnetic dipole pressure and 
using electron multiplicity (57), the maximum temperature is 

^^( 1 + |)" /2 - 9xl » 2 ( 1 + t)" /2 > < 58 > 

which applies as long as k&T p < m p c 2 . 

We therefore set the electron density to the threshold value for significant ion heating. Keeping 
in mind the results of §4.3, n e is fixed at a constant value (56) 

n e (Rmax) ( mpC 2^ 2Trm p c 2 (iZmag) *~ 59 ' ) 

along each field line (labeled by the maximum radius R mSLX ). The corresponding multiplicity pa- 
rameter * is given by eq. (57) or, equivalently, by 

A mag = ^(iU) (^) 3 - (60) 

The factor of (-R m ax/-Rmag)~ 3 m this expression accounts for the normalization of the multiplicity 
Amag at the magnetospheric boundary -R mag . 

One also expects a very strong dependence of T e on r when the electrons absorb the energy 
of large-scale torsional motions in the magnetosphere, and cool by incoherent synchrotron emission 
(SS4.2-4). The equilibrium temperature of the absorbing electrons is expected to increase rapidly 
with radius (eq. [30]). The synchrotron optical depth given by eqs. (39) and (50), with electron 
multiplicity (57), drops precipitously with distance from pulsar B, 

r u (r) oc rn e B 2 /*(vT e y^ oc e A A/ \~^ r~ 97 / 9 . (61) 
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(Here ea is the energy density in torsional motions at radius r in the pulsar B magnetosphere, 
relative to the magnetic pressure at the magnetopause boundary.) One finds, generically, a sharp 
transition from large to small optical depth, because of the much higher heating rates in the outer 
magnetosphere. Taking ea to be constant gives the frequency scaling 

A(f>(v) ~ v- 15/97 . (62) 

We therefore choose a temperature profile T e (r) oc r 7 . After setting the magnetospheric radius 
to -R mag = 4 x 10 9 cm (Mb = 1 in eq. [4]), the only free parameter in this model is the normalization 
of T e , which is set to k^,T e = Wm e c 2 at r = 1.4 x 10 9 cm. The resulting eclipse light curves are 
displayed in Fig. 9. One obtains an equally good fit to the central parts of the eclipse, and a 
smoother fall and rise in the radio flux at ingress and egress. The eclipse duration has a weak 
frequency dependence (Fig. 10), and the depth of the eclipse is less dependent on frequency than 
in the constant-n e model (compare Fig. 6). 

A rapid increase in the temperature toward the edge of the magnetosphere (combined with a 
drop in n e to maintain pressure balance) can therefore provide an explanation for the duration of 
the observed eclipses. Indeed, the dependence of T e on r could be even stronger than we assume 
if the amplitude ea of the turbulence increases outward in the magnetosphere of pulsar B, and if 
there is enhanced plasma loss from the outer magnetosphere. 

It should also be noted that our model of electron heating involves the absorption of low- 
frequency (1 — 10 MHz) radio waves (§4.2). This means that r v will maintain the scaling (61), and 
continue to drop off rapidly with radius, even when the plasma is optically thin at GHz frequencies. 

7. Polarization 

We now consider polarization effects associated with the propagation and absorption of the two 
electromagnetic modes in the eclipse region. In a magnetized, relativistic plasma these modes are 
linearly polarized below the synchrotron peak frequency, v < fp C ak = iy B,e{k^,T e /m e c 2 ) 2 , even if the 
plasma is composed of electrons and ions (e.g. Sagiv et al. 2004). The linear absorption coefficients 
of the two modes are given by eq. (39). We present a sample calculation of the polarization of 
the transmitted radio pulses in our favored eclipse model, under the assumption that the incident 
radiation is unpolarized. In this case, the radiation can become polarized as it propagates through 
an absorbing medium. 

In an inhomogeneous plasma (where the density and/or the magnetic field are functions of 
position), it is more convenient to describe radiation transfer in terms of modes polarized along 
fixed directions in space. We will choose these reference polarization states to lie along the normal 
to the orbital plane (a), and within the orbital plane (b). The radiation field is then characterized 
by four Stokes parameters I, Q, U and V. The first two parameters can be re-expressed in terms 



-23- 



of the intensities in the two reference polarization states, / 
polarization fraction is II = \J\J 2 + Q 2 + V 2 /I. 



I a + I b and Q = I a - I b . The 



Our calculation of the transmitted polarization neglects the effects of synchrotron re-emission 
and Faraday rotation. The first effect can be safely neglected if one is interested in the pulsed 
component. Faraday rotation is absent in a pair plasma, while for electron-ion plasma at low fre- 
quencies it is suppressed by a factor ~ ^/^ pC ak compared with the effects of synchrotron absorption. 
One consequence of this is that the effects of mode tracking are also negligible. This means that 
the instantaneous polarization angle, x = \ tan^ 1 (C//Q), is not able to adjust adiabatically to the 
changing direction of the magnetic field (e.g. Thompson et al. 1994). 

Marginally, effects of limiting polarization may lead to generation of small circular polarization, 
which depends on the difference n 2 — n\ between the mode indices of refraction depends on the 
angle between k and B. In a relativistic thermal plasma n 2 — n\ depends on the angle between k 
and B; approximately 



W r(i/3) 



0.54 a,, 



2\ 2 



< "B,e (kBT e /m e c 2 ) 



(63) 



where a„ is the polarization-averaged synchrotron absorption coefficient (e.g. Sagiv et al. 2004). 
In this situation, the most important consequence of this phase shift is to generate finite circular 
polarization. The circular polarization parameter V turns out to be small, but we include it in the 
calculation for completeness. 



Within the above approximations, the equations of polarization transfer read 
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(e.g. Pacholczyk & Swihart 1970). Here, xb is the angle between the reference direction a and the 
projected magnetic field (at a certain point x along the line of sight). 

The Stokes parameters are plotted in Fig. 11 for our best fit eclipse model. The transmitted 
radiation is predicted to be strongly polarized in the deepest part of the eclipse, with a polarization 
fraction reaching 25%. (This is, unfortuntely, the time when the signal to noise ratio is smallest.) 
The polarization angle evolves smoothly throughout eclipse. One finds U < Q when II is largest, 
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which means that the transmitted polarization tends to be aligned with one of the two reference 
polarization directions. The main contribution to U and Q comes from regions of the magnetosphere 
which have finite transparency. 

We find that the net value of V is much smaller than that of U and Q, by one or two orders 
of magnitude, in the case where the background source is unpolarized. This is partly the result of 
a cancellation between the ingoing and outgoing parts of the ray trajectory. 

Polarization provides an independent test of the eclipse mechanism and the geometry of the 
system. The radio pulses of pulsar A are, in fact, strongly polarized (up to 50%: Ransom et 
al. 2004). In order to predict correctly the transmitted polarization, one needs to know the 
direction of polarization with respect to the orbital plane. This can be achieved using calibration of 
polarization angle, i.e. by referencing to a source with known position angle, and a determination 
of the orientation of the orbit on the sky (by combining the scintillation measurements with the 
anticipated proper motion measurements). 

8. Predictions and Comparison with Other Models 

8.1. Predictions of the model 

We now summarize the predictions of our model, and how further observations can be used to 
probe the geometry of the system and the properties of the eclipsing material. 

1. The impact parameter zq between the line of sight and the orbital plane is predicted to be 
\zq\ = 7.5 ± 0.7 x 10 8 cm. The corresponding orbital inclination is 90.55° ± 0.05° (or possibly 
89.45° ± 0.05°). This prediction is intermediate between the measurements of Cole et al. 
(2004) and Ransom et al. (2004). 

2. The spin fls of pulsar B is expected to undergo geodetic precession on a ~ 75 year timescale 
(Burgay et al. 2003). We have provided predictions for how the eclipse light curve will vary 
as a result. In particular, the orbital phase at which the radio flux reaches a minimum will 
shift back toward superior conjunction (Fig. 13). Since 4>q is not well constrained, a time for 
eclipse to become symmetric is between ~ 12 years (if \<j)Q + 90° | ~ 30° and Qb is drifting 
away from 4>q = —90°) and ~ 25 years (if 4>q is drifting toward —90°). 

3. High time-resolution observations are a sensitive probe of the distribution of plasma properties 
(density and temperature) on the closed field lines. If plasma is depleted from the outermost 
field lines, at high temporal resolution the flux should return to unity. On the other hand, 
if the absorbing plasma does not have a sharp truncation in radius, then the flux will not 
return to unity in all of the transparent windows. 

4. The eclipses must regain a strong frequency dependence at sufficiently high frequencies. The 
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critical frequency above which significant transmission occurs can be used to place tight 
constraints on the plasma density. The electron cyclotron frequency at a distance r ~ zq ~ 
7.5 x 10 8 cm from pulsar B is estimated to be VB,e ~ 3 GHz. One therefore expects the 
eclipses to develop a significant frequency dependence at higher frequencies. 

5. There are a number of definite predictions for the polarization of the transmitted radiation as 
discussed in §7. The propagation of the radio waves through the intervening magnetosphere 
is in the quasi-transverse regime, where the polarization eigenmodes are linear and Faraday 
rotation can be neglected. The rotation measure should therefore not change significantly 
during the eclipse, and the variation in the dispersion measure is predicted to be small, of 
the order of 10 13 cm -2 . 



Finally, we note that the model predicts some radio emission from the magnetosphere of pulsar 
B, but its precise amplitude and spectrum cannot be easily determined from first principles. A 
straightforward upper bound to the radio output is obtained for incoherent synchrotron radiation. 
The luminosity of electrons (and positrons) radiating at temperature T e can be no larger than the 
blackbody at this temperature: 

dL hh _ 2irk B T e v 3 2 22 3 / k B T e /m e c 2 \ f R mag \ 2 , 

This is too weak to be important, and so the emission must be of a coherent type such as we have 
described. 



8.2. Alternative Eclipse Models 

In this section we briefly contrast our favored eclipse model with other mechanisms. At first 
sight, an important clue is that the electron cyclotron frequency, measured along the line of sight, 
peaks at a value close to the observing frequency. Thus, it is tempting to associate the eclipses with 
absorption of the radio waves on non-relativistic particles. Indeed, we have outlined a mechanism 
for supplying a constant flux of non-relativistic electrons (and ions) on the closed magnetic field 
lines (in Appendix A). If the closed dipolar field lines of pulsar B are twisted through an angle 
Acf) ~ 1 radian, then the optical depth through the cyclotron resonance (of either type of charge) 
is T cyc ~ P^ 1 . Here /3\\ is the drift speed of the charges, in units of the speed of light. Thus, a 
sub-relativistic drift of the charges (e.g. with a speed [3\\ ~ VLr/c) would guarantee a significant 
optical depth at the cyclotron resonance. 

There are two difficulties with this very simple model. First, the duration of the eclipse has 
a weak frequency dependence. The eclipses extend out to a distance R± ~ ±1.5 x 10 9 cm from 
superior conjunction (in a direction parallel to the orbital plane). For example, if the magnetic 
moment of Pulsar B were directed toward the observer, then the maximum electron cyclotron 
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frequency along the line of sight would be 

v B ,e = 400 + lj MHz (66) 

at eclipse ingress (or egress). Here zq is the offset of the line of sight from the orbital plane. The 
inclination of the orbital plane with respect to the line of sight is ^ — i, where i is the inclination 
with respect to the plane of the sky. One then has 

z = (i - |) D AB (67) 

The corresponding expression is 

/ 2 \ - 2 /a 2 \ !/ 2 

u B , e = 400 A^ /4 ^ + lj (i^ 1 ) MHz (68) 

when the magnetic moment is oriented perpendicular to the orbital plane. Inverting these ex- 
pressions, one obtains R± as a function of vs,e for a fixed impact parameter. In both cases, it 

— 1/3 

is a stronger function of frequency than R± ~ v B e when 6/0. In addition, one observes that 
absorption must occur at a harmonic > 1 unless the torque parameter Mb > 1 (see eq. 4). 

Cyclotron absorption inside the magnetosphere of pulsar B has some advantages: the cross 
section is larger than it is for synchrotron absorption by relativistic particles, and the required 
particle density is thereby reduced. To have an appreciable optical depth at a radius r (and corre- 
sponding frequency VB,e{f)) the particle density should exceed the local Goldreich- Julian density 
by a factor (f^r/c) -1 ~ 10. Second, the eclipse is observed to begin roughly where the line of sight 
passes deep enough into the magnetosphere of pulsar B that fB,e ~ v. Cyclotron absorption is, 
nonetheless, strongly disfavored for at least two reasons: first the expected frequency dependence 
of the eclipse duration is A<p oc z/ -1 / 3 (or a stronger function of frequency), in contradiction with 
the much weaker observed scaling Ac/) oc i/~ 0A ; and, second, the absorbing particles will inevitably 
be heated to relativistic temperatures. In addition, a fluctuating component of the current would 
cause a heavy fossil disk of cold particles, which could be centrifugally suspended in the outer 
magnetosphere, to be drained on a very short timescale (cf. Thompson et al. 2002). 

Our favored plasma parameters differ substantially from those advocated by Rafikov &: Gol- 
dreich (2004), basically because we invoke a much stronger heating mechanism, and a different 
source of relativistic particles. At high optical depth, particles are isotropized much more easily, 
and the equilibrium particle density regulates to a much larger value (due to the competition be- 
tween injection from below and precipitation through the cyclotron cooling radius). If the external 
radio photons are the primary heat source, as advocated by Rafikov & Goldreich (2004), then the 
optical depth at the cyclotron fundamental is regulated to value of the order of unity. One infers 
that the eclipse duration should vary as A</> ~ i^ -1 / 3 (or as a strong function of u if the radio pump 
radiation is brighter at ~ 100 MHz than it is at ~ 1 GHz). 
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We have also considered effects of induced Compton scattering of pulsar A radio emission by 
plasma in pulsar B magnetosphere and by pulsar B radio beam (e.g. Thompson et al. 1994). Both 
effects are negligible due to a small optical depth. 

9. Conclusion 

We have developed a model of the radio eclipses of pulsar J0737— 3039A, in which synchrotron 
absorption on relativistic particles occurs inside the intervening magnetosphere of the companion 
pulsar B. We believe that the value of such modeling is threefold. First, it provides a strong test 
of the longstanding assumption that isolated neutron stars are surrounded by corotating, dipolar 
magnetic fields. Second, one is probing how a pulsar magnetosphere interacts with an external wind, 
and in particular the mechanism by which turbulence is damped in a relativistic and magnetically 
dominated plasma. Third, an understanding of how charges are supplied to the magnetosphere 
may have broader implications for the electrodynamics of radio pulsars. 

Given the simplicity of the model, it is in excellent agreement with observations, especially in 
the middle of the eclipse where the interaction with the external wind has only a modest effect on 
the shape of the poloidal field lines. The model can explain the asymmetry of the eclipse between 
ingress and egress; the weak frequency dependence of its duration; the modulation of the pulsar 
A emission at both the spin period of pulsar B, and half its spin period, in different parts of the 
eclipse; the detailed shape of the luminosity spikes in the middle of the eclipse; and the dependence 
of the eclipse duration on the rotational phase of pulsar B. The model implies that pulsar B is 
nearly, but not exactly, an orthogonal rotator. It also requires the spin axis of pulsar B to be tilted 
from the normal to the orbital plane. 

We have demonstrated that - at intermediate distances - the poloidal magnetic field of a 
neutron star is well approximated by a dipole. This is a valuable confirmation of a fundamental 
assumption made in models of pulsar electrodynamics. In addition, our results are consistent with 
models that place the source of the radio emission close to the magnetic axis. Note that McLaughlin 
et al. (2004) define phases with respect to the arrival of radio pulses from pulsar B, whereas we 
define them with respect to the magnetic axis of B. On the other hand, the dependence of eclipse 
duration on pulsar B phase in our model is the same as inferred by McLaughlin et al. (2004), so 
that the two definitions of phases are close to each other. 

In the present model we have neglected the fact that the radio pulse profile of pulsar B has a 
single peak. This can be used to put addition constraints on the geometry of the system, assuming 
that the radio emission mechanism is the same at both magnetic poles. In particular, in order for 
the radio beam from one of the pole to miss the observer, the angle (f)Q should not be equal to 
±90°. This comes from the fact that for c/)q = ±90° the angle between the direction of the magnetic 
moment and the line of sight is the same every half a period. Unfortunately, in the absence of pulsar 
radio emission model with predictive power, this does not provide any meaningful constraint at the 
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moment. 

Our eclipse model has implications for the formation and evolution of the PSR J0737— 3039A 
system. The tilt of the spin axis of pulsar B with respect to the normal to the orbital plane is 
consistent with a large natal kick (as has been inferred from the spatial velocity of the system: 
Ransom et al. (2004)). The kick would change the orientation of the orbital plane and disrupt 
any pre-existing alignment between orbit and the spin of the progenitor star. Our measurement 
also does not support the suggestion of Demorest et al. (2004) that the spin axis should become 
aligned with normal to orbital plane due to torque from pulsar A. 

We have also described how an optically thick, relativistic, thermal plasma may be maintained 
in the outer magnetosphere of pulsar B, where synchrotron cooling is slow. We have argued that 
the interaction between the wind of pulsar A and the magnetosphere is much more effective at 
heating the particles than is the absorption of the radio pulses of pulsar A. This interaction also 
generates two sources of particles on the closed magnetic field lines. First, bunches of relativistic 
charges are created by a pulsar-type mechanism close to the neutron star surface as the dipolar field 
lines are twisted back and forth close to the magnetic poles. Second, magnetic helicity builds up on 
the dipolar field lines, as the result of the asymmetry between the outgoing and return currents in 
each magnetic hemisphere. The resulting steady current supplies a constant flow of electrons and 
ions from the star to the outer magnetosphere. 

The relativistic charges flowing in the magnetosphere of pulsar B are generally a more potent 
source of radio photons than is the companion pulsar. The requirement for this to be true is 
that the efficiency of conversion of turbulent energy to radio photons is larger than ~ 10 . (For 
example, if the unpulsed radio emission originates in the magnetosphere of pulsar B, then this 
minimal efficiency is exceed by ~ 10 4 .) As a result, heating of the magnetospheric particles occurs 
mainly through self-absorption of these internally generated photons. Heating of ions can have 
an important regulatory effect on the equilibrium electron density: if n e exceeds a critical value, 
then radio photons that resonate with the ion gyromotion are absorbed by the plasma. A plausible 
source of low frequency pump photons has also been identified. This involves a two-step emission 
process: a cascade of torsional Alfven waves to high frequencies where the waves become charge 
starved; followed by wiggler emission of radio waves by the electrostatically accelerated charges in 
the fluctuating magnetic field. 

We would like to thank Maura McLaughlin for sharing original data and comments on the 
manuscript. We also would like to thank Roman Rafikov, Peter Goldreich, Jeremy Heyl, Victoria 
Kaspi, Anatoly Spitkovsky, Ingrid Stairs, and Scott Ransom for stimulating discussions. 
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A. Helicity Pumping on the Closed Magnetic Field Lines 

The excitation of large-scale torsional motions provides a robust mechanism for drawing charges 
into the outer magnetosphere of pulsar B. Damping of this turbulence will also heat the suspended 
particles (§4.2), and if the heating rate is too high then the plasma may not absorb radio waves 
effectively. This leads us to examine a subtler effect, involving a static component of the magne- 
tospheric twist, which will operate deeper in the magnetosphere of pulsar B where the fluctuating 
current is small, \SJ\ <C |pgj|c 

The exterior magnetic field of a neutron star can store a considerable amount of magnetic 
helicity. When it does, the current flowing nearly parallel to the magnetic field has a zero-frequency 
component. In the magnetically active Soft Gamma Repeaters and Anomalous X-ray Pulsars, the 
unwinding of a strong, internal magnetic field is a repeating source of magnetic helicity for the 
exterior of the star (Thompson et al. 2002). One does not expect such a mechanism to be active 
in the much older (and weakly magnetic) PSR J0737— 3039. 

Nonetheless, we offer a simple argument why magnetic helicity will be pumped from the outer 
to the inner magnetosphere of pulsar B. There are two simple components to this argument. First, 
if a fixed amount of helicity is injected into a dipole magnetic field, then the energy in the toroidal 
magnetic field is minimized if the helicity 

H ~ {ir^fBpB^ ~ (TtB P r 2 ) 2 A(f) (Al) 

is concentrated close to the star (Thompson et al. 2002). Here Bp is the poloidal magnetic field 
at radius r, and B^ ~ BpA(f>. The energy in the toroidal magnetic field is 

e *~^'- 3 ~isSij kwV < A2) 

If a magnetosphere carrying helicity TL is not in the minimum energy state, the winding of the 
magnetic field can be rapidly transferred by reconnection between different flux surfaces (e.g. Taylor 
1974). 

The second observation is that the current flowing out along an open bundle of magnetic field 
lines will not be compensated precisely by the return current in the same magnetic hemisphere. 
There is, as a result, a residual current connecting the two magnetic poles, which flows along the 



This preprint was prepared with the AAS IATgjX macros v5.2. 



-31 - 



outermost closed magnetic field lines. This closed field current imparts a small amount of helicity 
to the magnetosphere. By the above argument, this helicity will be pumped, by reconnection, on 
a short timescale ~ c/f2 to the inner magnetosphere. The net rate of transfer of helicity is then 

dH 



dt 



e n $l pcn n. (A3) 



Here <I>open is the magnetic flux $(r) = irBp(r)r 2 evaluated at the boundary radius of the magne- 
tosphere. 

In the case of an isolated neutron star, an imbalance between the outgoing and return current 
can result from an offset AR^ po \ c of the magnetic dipole from the center of the star. One then 
deduces 

f2Ai? dipolc S7i? NS . . . 

£ H r V— < . (A4) 

c c 

Making use of 3> pen = ^-B P o1c-Rns(^^ns/ c), one can relate the helicity increase (A3) to the mag- 
netic dipole luminosity L M dr = ^{^R-ns/c) a B^ o1c R^ s c, 

~ 67r 2 L M DRAi?di P oic- (A5) 

+ 



dt 



where -B po i e is magnetic field at the magnetic pole of the star. 

The effect is much stronger in the case of PSR J0737— 3039B. The strong asymmetry in the 
external stress of the pulsar A wind will enforce a large asymmetry between the outgoing and return 
currents in each magnetic hemisphere. We therefore estimate 

e H ~ A^iCag) ~ ~ 0.1 - 0.2. (A6) 

The equilibrium twist angle that is established at smaller radii depends on how rapidly the 
helicity is damped in the inner magnetosphere. We will not examine this problem here, and simply 
note that the inward flux of helicity allows one to establish a minimal static twist angle Acfr in the 
outer parts of the magnetosphere. Given that the twist is transferred between flux surfaces with a 
speed ~ c, one has 

dH 



dt 

Combining this with eqs. (A3) and (A6) gives 



§\r) A<f>(r)-. (A7) 



r 



Aco(r)~('^^) (-— ) . (A8) 



The corresponding current is 

j f r y (A9) 

IpgjIc V 0.3 J \R ma J 
It will be noted that the presence of a zero-frequency component of the current will influence the 
non-linear couplings between torsional waves in the magnetosphere (Goldreich & Sridhar 1997; 
Thompson & Blaes 1998). 
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Fig. 1. — Geometry of the model. The absorbing plasma is axisymmetrically distributed about 
pulsar B, and extends as far as a distance i? a bs+ < -Rmag from the neutron star. The spin axis 
of pulsar B is defined by angles 9q and <j)Q. The magnetic moment /7b makes angle x with the 
spin axis and rotates around tts every 2.77 s, as the pulsar is carried by the orbital motion in the 
y-direction. The line of sight passes at distance zq above the orbital plane. 
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Fig. 2. — View of the magnetosphere at different rotational phases separated by 7r/4. For a full 
movie of the eclipse see http://www.physics.mcgill.ca/~lyutikov/movie.gif 
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Fig. 3. — Inferred geometry of the pulsar B magnetosphere projected onto the z — y plane at 
the two particular moments where jls lies in the z — y plane. The spin axis fig is inclined by 
angle ~ 60° to the z axis and lies in the y — z plane. The magnetic moment \ib makes an angle 
~ 75° with fts- Apparent velocity of pulsar is along —y axis. The line of sight passes at 
zq/ 'i? a bs+ ~ —0.5. The absorbing plasma is located between the two drawn magnetic surfaces with 
i? a bs- = -75i? a bs+ < R < -R a bs+- The inner circle marks the cooling radius -R C ooi = -05i2 a b s +- The 
dotted line is the equatorial plane, orthogonal to fie. The two directions of hb plotted correspond 
to pulsar B phases separated by half a period. Points corresponding to negative values of y enter 
the line of sight first. The transition to one transparent window per rotation occurs near the point 
when the line of sight crosses the dashed field lines. A full eclipse occurs when the line of sight 
crosses both solid and dashed field lines. 
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Fig. 4. — Comparison of a simulated eclipse profile (solid line) with 800 MHz data (McLaughlin 
et al. 2004). In the model, the absorbing electrons are distributed on the dipole magnetic field 
lines with a constant density A ma g = 3 x 10 5 (/c B T e /10m e c 2 ) 5 / 3 ^Hz over a fixec ^ ran S e of ra dh 
Rabs- < flmax < -R a bs+- Here i? a bs+ = 1-5 x 10 9 cm and i? a bs- = 0.75-R abs +- The same light 
curve will result if v, T e and A mag varied jointly so as to keep A mag /(T e i/) 5//3 fixed. The spin axis 
of pulsar B is oriented at Oq/(^tt) = 0.66 and 4>q = —\tt. The inclination between Qb and the 
magnetic moment /7b is x/{\ 7r ) = 0.82, and the impact parameter is zq = —7.5 x 10 8 cm. The 
model light curve is smoothed over a timescale At = 0.25 s with a gaussian window e~ A ^ ^ 2<T « > 
[cj^ = At/(27r) 1 / 2 ]. The model fits the data best in the middle of the eclipse. 
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Fig. 5. — Same as Fig. 4, but model and data side by side. 
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Fig. 6. — Frequency dependence of the eclipse profile, smoothed on a timescale At = 2.77 s. 
Parameters are the same as in Fig. 4. At high frequencies eclipse duration should become strongly 
frequency-dependent. This can be used to place tighter constraints on plasma density. 
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Fig. 7. — Average eclipse profiles for different orientations of magnetic moment of B. Solid line: 
eclipse profile averaged over points 315° < 4>b < 45° and 135° < (ps < 225°, where 4>b is a 
rotational phase of B (efts ~ 0°, 180° corresponds to magnetic moment pointing nearly towards an 
obsever, while 4>b ~ 90°, 270° corresponds to magnetic moment in the plane of the sky). Dotted 
line: eclipse profile averaged over 45° < 4>b < 135°, dashed line: eclipse profile averaged over 
225° < 4>b < 315°. Eclipse is broader when the pulsar is "face-on". Parameters are the same as in 
Fig. 4. 
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Fig. 8. — Change in the eclipse light curve due to a variation in the impact parameter zq and the 
inclination \ between the magnetic and spin axes of pulsar B. Parameters are as in Fig. 4. 
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Fig. 9. — Comparison of a simulated eclipse profile with 800 MHz data (McLaughlin et al. 2004) 
in the model with large ion content (§6). The density depends on radius according to Eq. (57), 
and the temperature distribution is k^,T e = 10 m e c 2 (r/l. 4 x 10 9 cm) 7 . The magnetospheric 
radius is -R mag = 4 x 10 9 cm (Mb = 1 in eq. [4]). The other parameters z®, Qb, and \ are the same 
as in Fig. 5. The model light curve is smoothed over At = 0.25 s. 
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Fig. 10. — Frequency dependence of the eclipse profile in the model of Fig. 9, smoothed on a 
timescale At = 2.77 s. Parameters are the same as in Fig. 9 
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Fig. 11. — Polarization generated during the eclipse of an unpolarized background source, in the 
model of Fig. 4. The reference directions for the linear polarization are the normal to the orbital 
plane (a); and the orbital plane (b). Only negligible circular polarization is produced (V ~ 0). 
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Fig. 12. — (a) Predicted eclipse profile at high temporal resolution of At = 0.01 s at 800 MHz. 
Parameters as in Fig. 4. The absorbing plasma has a sharp cutoff on dipole field lines extending 
beyond -R a bs+ = 1-5 x 10 9 . In the transparent windows, the flux returns to the unabsorbed level, 
(b) Eclipse model with electron density given by eq. (56). Parameters as in Fig. 9. In some 
transparent windows, the flux is partially absorbed. 
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Fig. 13. — Prediction of the effects of geodetic precession on the eclipse light curve. The parameters 
are the same as in Fig. 9, but with a temporal resolution of At = 0.01 sec. 



